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Abstract. We study the homogenization of a diffusion process which takes 
place in a binary structure formed by an ambiental connected phase surrounding 
a suspension of very small spheres distributed in an e-periodic network. The 
asymptotic distribution of the concentration is determined for both phases, as 
e — > 0, assuming that the suspension has mass of unity order and vanishing 
volume. Three cases are distinguished according to the values of a certain limit 
capacity. When it is positive and finite, the macroscopic system involves a 
two-concentration system, coupled through a term accounting for the non local 
effects. In the other two cases, where the capacity is either infinite or going to 
zero, although the form of the system is much simpler, some peculiar effects still 
account for the presence of the suspension. 
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1 Introduction 



Diffusion occurs naturally and is important in many industrial and geophys- 
ical problems, particularly in oil recovery, earth pollution, phase transition, 
chemical and nuclear processes. When one comes to a rational study of binary 
structures, a crucial point lies in the interaction between the microscopic and 
macroscopic levels and particularly the way the former influences the latter. 
Once the distribution is assumed to be e-periodic, this kind of study can be 
accomplished by the homogenization theory. 

The present study reveals the basic mechanism which governs diffusion in 

both phases of such a binary structure, formed by an ambiental connected phase 

surrounding a periodical suspension of small particles. For simplicity, the par- 

v 

tides are considered here to be spheres of radius r e << e, that is lim — = 0. 

£->0 £ 

We balance this assumption, which obviously means that the suspension has 
vanishing volume, by imposing the total mass of the suspension to be always 
of unity order. This simplified structure permits the accurate establishment of 
the macroscopic equations by means of a multiple scale method of the homoge- 
nization theory adapted for fine-scale substructures. It allows to have a general 
view on the specific macroscopic effects which arise in every possible case. As we 
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use the non-dimensional framework, the discussion is made in fact with respect 
to only two parameters: r e and b e , the latter standing for the ratio of suspen- 
sion/ambiental phase diffusivities. As the diffusivities of the two components 
can differ by orders of magnitude, the interfacial conditions play an important 
role. 

It happens that the following cases have different treatments: r £ << e 3 , 
e 3 << r e << e and r e = 0(e 3 ). 

To give a flavor of what may be considered as an appropriate choice of the 
relative scales, we refer to the pioneering work [7J where the appearance of an 
extra term in the limit procedure is responsible for a change in the nature of the 
mathematical problem and is linked to a critical size of the inclusions. Later 6 
showed how this could be generalized to the iV-dimensional case for non linear 
operators satisfying classical properties of polynomial growth and coercivity. 
Since then, the notion of non local effects has been developed in a way that is 
closer to the present point of view in @|, and [HJ- 

In dealing with our problem, the main difficulty was due to the choice of test 
functions to be used in the associated variational formulation and which are 
classically some perturbation of the solution to the so-called cellular problem. 
Indeed, proceeding as usual in homogenization theory, we use energy arguments 
based on a priori estimates where direct limiting procedure apparently leads 
to singular behavior. Non local effects appear when these singularities can be 
overcome, which is usually achieved by using adequate test functions in the 
variational formulation. Since the fundamental work [7J, an important step 
was accomplished in this direction in 0. A slightly different approach [H] uses 
Dirichlct forms involving non classical measures in the spirit of ^H]- However, 
the main drawback of this method lies in its essential use of the Maximum Prin- 
ciple, which was avoided in 0] for elastic fibers, and later in [5] where the case of 
spherical symmetry is solved. The asymptotic behavior of highly heterogeneous 
media has also been considered in the framework of homogenization when the 
coefficient of one component is vanishing and both components have volumes of 
unity order: see the derivation of a double porosity model for a single phase flow 
by [2] and the application of two-scale convergence in order to model diffusion 
processes in £Q. 

The paper is organized as follows. Sectionals devoted to the main notations 
and to the description of the initial problem. We set the functional framework 
lfTH|l where the existence and uniqueness of the solution can be established: see 
|H] and for similar problems. In Section we introduce specific tools to 
handle the limiting process. This is based on the use of the operators G r defined 
by H38fl which have a localizing effect: this observation motivates the additional 
assumption (|46|l on the external sources when the radius of the particles is of 
critical order e 3 with e denoting the period of the distribution. While passing to 
the limit, the capacity number 7 £ defined by (|33fl appears as the main criterium 
to describe the limit problem, the relative parts played by the radius of the 
particles and by the period of the network becoming explicit. 

Section which is actually the most involving one, deals with the critical 
case when 7 e has a positive and finite limit 7. In this part, where we assume 
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also b £ —>■ +00 the test functions are a convex combination of the elementary 
solution of the Laplacian and its transformed by the operator G Te defined by 
(|38[) with r — r £ . This choice, which is inspired from [3], 0] and [8] and has to be 
compared with 7 , allows to overcome the singular behavior of the energy term 
when the period e tends to zero. We have to emphasize that this construction 
highly depends on the geometry of the problem, that is the spherical symmetry 
. To our knowledge, the generalization to more intricate geometries remains to 
be done. The resulting model (|69|l -l|72 () . with the initial value defined after uo 
in l|21[l and vq in H23|) . involves a pairing (u, v) which is coupled through a linear 
operator acting on the difference u — v by the factor 4ir"f. 

The case of the infinite capacity, where e 3 << r e << e, is worked out in 
Section The proofs are only sketched because the arguments follow the same 
lines as in Section 0] Let us mention that the singular behavior of the capacity 
in this case, that is j e — * +00, forces v to coincide with u. In other words, 
the infinite capacity prevents the splitting of the distribution, as it did in the 
critical case. Quite interestingly, the initial value of the global concentration 
is a convex combination (|83[) of the initial conditions uq and vq; moreover, the 
mass density of the macroscopic diffusion equation l|82[) takes both components 
into account, in accordance with the intuition that the limiting process must 
lead to a binary mixture. 

Finally, the case of vanishing capacity is handled in Section [SI that is when 
r e << e 3 . Here, v remains constant in time, obviously equal to the initial condi- 
tion vq, while u satisfies the diffusion equation JHJJ-JSSJ with data independent 
of the initial condition of the suspension. This can be seen as a proof that when 
the radius of the particles is too small, then the suspension does not present 
macroscopic effects, although a corresponding residual concentration, constant 
in time, should be considered. 



2 The diffusion problem 



We consider ft C R 3 a bounded Lipschitz domain occupied by a mixture of two 
different materials, one of them forming the ambiental connected phase and the 
other being concentrated in a periodical suspension of small spherical particles. 
Let us denote 

Y £ k :=ek + eY, keZ 3 . (2) 
Z £ := {k e Z 3 , Y £ k C n}, n Yc := U keZ Y £ k . (3) 
The suspension is defined by the following reunion 

D e :=U keZs B(ek,r e ), (4) 

where < r e << e and B(ek, r s ) is the ball of radius r e centered at ek, k £ Z e . 
Obviously, 

|L» e |^0 as e^O. (5) 
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The fluid domain is given by 

n e = n\D E . (6) 

We also use the following notation for the cylindrical time-domain: 

n T := fix]0,T[; (7) 

similar definitions for SlJ , Sly and Dj. 

We consider the problem which governs the diffusion process throughout our 
binary mixture. Denoting by a £ > and b e > the relative mass density 
and diffusivity of the suspension, then, assuming without loss of generality that 
|f2| = 1, its non-dimensional form is the following: 

To find u £ solution of 

p £ ^~dw{k £ Vu £ ) = f in n T (8) 

[u% = on dD T E (9) 

[k £ Vu%n = on dD T E (10) 

u £ = on dn T (11) 

u £ (0) = u £ in fi (12) 

where [-] e is the jump across the interface dD £ , n is the normal on dD e in the 
outward direction, f £ e L 2 (0, T; e L 2 (f7) and 

'«-{!. S«2; < 13 > 

Let H £ be the Hilbert space L 2 (Q) endowed with the scalar product 

{u,v) He := (p £ u,v) n (15) 
As Hq(D.) is dense in H e for any fixed e > 0, we can set 

#o( n ) QH e ~H' e Cif- x (0) (16) 

with continuous embeddings. 

Now, we can present the variational formulation of the problem <|8|)- (|12|) . 

To find u £ e L 2 (0,T;H£(tt)) n L°°(0, T; ff £ ) satisfying (in some sense) the 
initial condition 1|12[) and the following equation 

^(u e ,w) He +(fc £ Vu £ ,V W )o = (f e ,w) in 2/(0, T), Vu> G -ffg(il) (17) 
where (•,•) denotes the duality product between ff _1 (f2) and Hq(Q). 
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Theorem 2.1 Under the above hypotheses and notations, problem ji7| ) has a 

unique solution. Moreover, — — £ L 2 (0,T; H^ 1 ^)) and hence, u e is equal 

almost everywhere to a function ofC ([0,T];H e ); this is the sense of the initial 
condition \12]) . 

In the following we consider that the density of the spherical particles is 
much higher than that of the surrounding phase. The specific feature of our 
mixture, which describes the fact that although the volume of the suspension is 
vanishing its mass is of unity order, is given by: 

lim eUlU = a > (18) 
Regarding the relative diffusivity, we only assume: 

b £ >b Q > 0, Ve > 0. (19) 

As for the data, we assume that there exist / £ L 2 (0, T; iJ _1 (J7)) and uo £ 
L 2 (fl) such that 

f^f in L 2 (0,T;H-\n)) (20) 

u e u in L 2 {tt) (21) 

Also, we assume that there exist C > (independent of e) and vq £ L 2 (U) for 
which 

\u e \ 2 dx < C (22) 



' u e XD s v in V\n) (23) 



1^. 

where, for any D C 11, we denote 

/• dx = -T-— / -dx. 
d \D\ J D 

Remark 2.2 As Uq satisfies then holds at least on some subsequence 
(see Lemma A-2 

Proposition 2.3 We have 

u e is bounded in L°°(0,T; L 2 (Q)) n L 2 (0,T; H^(fl)). (24) 

Moreover, there exists C > 0, independent of e, such that 



e|2 

I'" 

D, 



-dx < C a.e. in [0,T] (25) 



be\Vu e \ 2 L2(DI) < C. (26) 
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Proof. Substituting w — u e in the variational problem l|17l) and integrating 
over (0, t) for any t G]0, T[, we get: 



t rt 



\(\u e (t)\^+a e \u e {t)\l s )+b e I \Vu%ds+ \ \Vu%Js 



JO 



= I (/ e (*),« e («)>ds + i(|«§|^+a e |uS|?,.) 



Notice that |2TJ and © yield: 

l«oln. + *eWa\l e < Kin + ^\D e \ I \u%\ 2 dx < C. 

J D e 



Moreover: 



/ {f e {s),u e {s))d S < f \f\ H -^\nd S 
J o Jo 

< I \f\H-i\Vu 6 \n e ds+ f \f\ H -r\Vu e \ De ds 
Jo Jo 

< \ J]f%-ids + \ f\Vu%ds + ±- jjril-.ds + | f\Vu% E ds. 

There results: 

~ (|« e (t)|^ + a e \u%t)\ 2 D J + | J* \Vu% E ds + \ j* \Vu%ds < C 
and the proof is completed. ■ 

3 Specific tools 

First, we introduce 

K £ = {R, r £ « R « e} 

that is R £ TZ £ iff 

lim — = Jim - = 0. (27) 
e^o R e^o e K ' 

We have to remark that lZ e is an infinite set, this property being insured by 
the assumption < r £ << s. 

We denote the domain confined between the spheres of radius a and b by 

C(a,b) := {x e R 3 , a < \x\ < b} 

and correspondingly 

C k {a,b) :=efc + C(a,6). 
For any R £ £ 7Z £ , we use the following notations: 

C £ :=U ke z E C k (r £ ,R £ ), C e T := C E x]0, T[ 
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Definition 3.1 For any R £ E 1Z £ , we define w Rc E Hq(CI) by 



in il e \C e , 
w Rs (x) := { W Re (x-sk) in C e fe , Vfc E Z e , 

1 in D F . 



where 



w Re (y) 



R, 



1 for yEC{r e ,R e ) 



(28) 



(29) 



{Re ~ r e ) V Itfl 

We have to remark here that W Re E H {C(r e , R £ )) and satisfies the system 

AW Rc = in C{r E ,R E ) 
W Rs = 1 for \y\ = r e 
W Rc = for \y\=Re 

From now on, we denote 

Proposition 3.2 For any R e E TZ S , we have 

|Vu^|n<C 7e 1/2 
w Rs — > in L 2 (Sl). 

Proof. First notice that 



(30) 
(31) 
(32) 

(33) 

(34) 
(35) 



\wr. \n = \w Rc \ Cc uD e < \C e UD e \^ 2 <C 

and lim e ^o — = by assumption (|27|) . 
As for the rest, direct computation shows 

\V«>R.\h = E / |Vwfl.| 2 dx 

= E 



3/2 



C B n_ i 



J<I> / sin6d6 
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7 \R E -r E 



< C 



Re-r e 

le 



and the proof is completed by (|27|l . ■ 

Lemmas 13.31 and 13.41 below are set without proof since they are a three- 
dimensional adaptation of Lemmas A. 3 and A. 4 0. 

Lemma 3.3 For every < r\ < r2 anc! u E H[C{r 1^2)), the following esti- 
mate holds true: 



C(n,r 2 ) 



> 



47rrir2 
^2 - ri 



+• u da — 


j- u da 


J Sr 2 


1 s ri 



(36) 
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where 

1 



da := - — - / • da 



s 4?rr 2 Js 



Lemma 3.4 There exists a positive constant C > such that: V(R, a) S R + x 
(0,1), Vu G H x (B(0,R)), 

[ \u~ I uda\ 2 dx<C—\Vu\ 2 B(aR) . (37) 

JB(0,R) J S aR OL 



Definition 3.5 Consider the piecewise constant functions G r : L 2 (0, T; (il)) 
L 2 (Q T ) defined for any r > by 

G r {6){x,t) = V ( / lr*W (38) 

feez e V s* 7 

where we denote 

S k r =dB(ek,r). (39) 



Lemma 3.6 If R e E TZ £ , then for every 6 <E L 2 (0, T; Hq (fi)) we Ziaue 

-3 \ V2 



L 2 (n T ) 



(40) 



|fl-G re (0)| iW) < Cr E |V0U W) (41) 

~3\ 1/2 

|V0| L 2 (C T) (42) 



\GrM-g t M\l^) < c(^j 

where GR e (6) and G re {0) are defined following f .¥<**)) . 
Moreover: 

\GrM\1^) =£-fjG R M\ 2 dxdt, \G r M\l H nT } = £f JG re (e)\ 2 dxdt. 



(43) 



Proof. Notice that by definition: 

i-T r r rT 



I I \ 9 ~{ 9da\ 2 dxdt<J" f f \0--f 6 da\ 2 dxdt 

C7l Jo Jy* J si h J B {ek,^.) J s* 



where we have used that 

for every k g Z e . We use Lemma f3. 41 with 

if = — — , a 



V3 
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to deduce that 



/ \6- G Re (9)\ 2 dxdt < C I ^) £ [ T [ | W | 2 dxdt 

C^- [ [ _ \V9\ 2 dxdt < C^- I \V9\ 2 dxdt 

K-z Jo JBIek,^) K e Jqt 



< _ 

which shows H40(l . 

To establish ff 4 1 1> - we recall the definition: 

f \9-G rE {9)\ 2 dxdt = V / / |0-/ 9da\ 2 dxdt 

Applying Lemma 13 . 41 with R — r e and a = 1, we get the result 
/ \0- G r X0)\ 2 dxdt < Crl V / / |V6>| 2 dxdt < Cr 2 [ |V6»| 2 dxdt. 

JDT keZ ^Jo JB(ek,r c ) J DJ" 

We come to (021- Indeed, applying Lemma 0751 and i|27|): 

/ |G B (0)-G P (6>)| 2 dzdt = V // |/ 9 da - -f 9 da\ 2 dydt 

kez e J ° J n J s k He J 

= Ce 3 < ^ LZ rr^ I |V^| 2 dxdt < C— [ \V9\ 2 dxdt. 
Finally, a direct computation yields (|43|l . ■ 

Proposition 3.7 If R e £ K E , then for any 9 E L 2 (0, T; Hq(Q)) there holds 
true: 



f I \9\ 2 dxdt < Cmax(l,— )|V0| 

JO J D E r e 



3 

2 

L 2 (Q. T ) 



Proof. We have 

i-T r t-T 



[ I \9\ 2 dxdt<2[ I \e-Gr E {9)\ 2 dxdt + 2 [ I \G rE {9)\ 2 dxdt 

J0 J D E Jo J D € Jo J D e 

= 2 [ / \9-G rE (9)\ 2 dxdt + 2 [ \G rc {9)\ 2 dxdt 
Jo J d, Jn T 



<Cr 2 f / \V9\ 2 dxdt + Af \G rE {9) - G Rc {9)\ 2 dxdt+ 
Jo J D e Jn T 

+8 / \G Re (9) - 9\ 2 dxdt + 8 [ \9\ 2 dxdt 
Jn T Jn T 



< c 



<Cr\ [ f \V0\ 2 dxdt + C— [ |V6>| 2 dxdt+ 

Jo J D e r e Jc? 

+C— [ |V(9| 2 dxdt + C [ |V6»| 2 dxdt 
Re Jqt Jn T 

— + — + 1 / |V6»| 2 dxdt < Cmaxfl,— ) / \W9\ 2 dxdt 

r e Re J Jqt r £ J QT 



Remark 3.8 Using the Mean Value Theorem, we easily find that 

\GrAf) ~ <P\L"°(C s \JD e ) < 2i?e|V^| L oo (n) , V(fi £ T>(D,). 

Definition 3.9 Let M De : L 2 (0, T; C C (Q)) -> L 2 (n T ) be defined by 



M Da (<p){x,t) := J2 ( / 



<P(V,t) dy lB{ek,r s ){x). 



Lemma 3.10 For any <p G £ 2 (0, T; C c (fi)), we /lave: 



Proof. Notice that 



lini / / \ip- M De (tp)\ 2 dxdt = 0. 

' i D r 



/ ^-M^O^da^-^- V / |^-/ <pdy\ 2 dx. 

•/ -D E l-^el t,- JB(ek,r e ) J Y, k 



k£Z s J B{ek,r e ) 

ini , , ,,card(Z £ ) 

y, then |B( ,r e )|^- £ 

continuity of y> on £1 it follows the convergence to a.e. on [0, T]. Lebesgue's 
dominated convergence theorem achieves the result. ■ 



As card(Z £ ) ~ then \B(0,r e )\ — , n , -> |0| = 1 and by the uniform 



4 Homogenization of the case r £ = 0(s 3 ) 

The present critical radius case is described by 

lim7 e = 7 e]0, +oo [. (44) 

e^O 

Its homogenization process is the most involving one. That is why we start the 
homogenization study of our problem with this case, under the condition 

lim b £ = +oo (45) 

We also assume that f £ has the following additional property: 
3R E G TZ e and g G i 2 (0,T; for which 

(46) 

(f s ,w Re tp)^(g,<p) in Tf(0,T), V^6P(!1) 
(see |H] for a certain type of functions f £ which satisfy 14611 1. 
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Remark 4.1 Notice that due to \44V r i n this case Proposition \3. 7| reads 

VcpeL 2 (0,T;H^Sl)), f I \^\ 2 dxdt<C\S7^\ 2 L2{nT) . (47) 

JO J D c 

A preliminary result is the following: 

Proposition 4.2 There exist u E L°°(0, T; L 2 {£1)) n L 2 (0, T; Hq(Q)) and v € 
L 2 (Q T ) such that, on some subsequence, 

u £ ^u in L°°(0,T;L 2 (Q,)) (48) 

u £ ^u in L 2 {Q,T;Hl{n)) (49) 

G Rf {u £ )^u in L 2 (fl T ) (50) 

G re (u e )^v in L 2 (n T ) (51) 

Moreover, we have 

lim / / \u £ - G r {u £ )\ 2 dxdt = (52) 

Proof. From (I24|) . we get, on some subsequence, the convergences (|48|) and 
<|49|) . Moreover, we have: 

W-G Rc (u £ )\l T = \u\ 2 T +\ u -G RE (u £ )\ 2 nT (53) 

where: 

l«-G fle ( U £ )| n T <| U -w e |n- +K-G fle (?/)| T (54) 
< \u~u £ \ nT + \u £ ~G Re {u £ )\ n T 

and gOJl yields: 



and thus: 

£^0 



hmJ^-G^)!^ =0. 



As (03 implies that 

u £ ^u in L 2 (f7 T ) (55) 
the right-hand side of l|54|l tends to zero as e — > 0, that is: 

lim|u-G Re (<)| n? =0. 
After substitution into the right-hand side of (|53|l . and taking into account that 

lim \n T \ fi£ I = 0, 

we obtain JSUJ, that is, 

G Re {u £ )^u in i 2 (fi T ). (56) 
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In order to prove JHH, we see that 

|G r > E )U W < \G rc (u E ) ~ G Rc {u £ )\ LH nT) + \G Re (u e )\ LHQT) 

c (57) 
< —\S7u e \ v . ( ar ) + C<C. 

Moreover, recall that from (|41|) we have, taking into account f26|) : 



f I \u £ - G Te (u £ )\ 2 dxdt < Crl f f \Vu £ \ 2 dxdt < 

JO J D e JO J D e 

and the proof is completed. 



C 

7 A 



(58) 



Proposition 4.3 For any p £ L 2 (0,T;C c (fl)), we have: 

lim / + u £ ipdxdt = / vipdxdt. (59) 
Jo J D e Jn T 

Proof. We have: 

u £ ipdxdt= I j- (u £ — G rc (u £ ))ipdxdt+ 

(60) 



J D s JO J D 



+ 11 G rs {u £ ){ip~ M Ds (tp))dxdt+ [ I G rc {u £ )M Ds (ip)dxdt 

Jo J D e Jo J D e 

The first right-hand term tends to zero thanks to (|52|) in Proposition ^. 21 The 
second one tends also to zero thanks to Lemma \'S. 101 The last term is handled 
as follows: 

f I G rE (u £ )M DE (ip)dxdt = \ E V] ([ [I u £ da] if dxdt = A e / p>G Te {u £ )dxdt 
Jo J D e k£Z Jo Jy* \J J Jn T 

where 

The proof is completed by (|5T|l . ■ 
Proposition 4.4 For any p> G L 2 (0,T; Hq(Q)), we have 

u £ ipdxdt — > / vipdxdt. (61) 
Ja T 



Proof. In the light of proposition 14. 31 we have to prove that the left-hand side 
term is continuous in the corresponding norm. This can be obtained as follows: 

1/7- irVd.vdl I I I I U 2 dxdt\ < 



u £ pdxdt 
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I I \u £ \ 2 dxdt] ( I I \>p\ 2 dxdt] 

^Jo J D E J \Jo J D e J 
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< C|^lL 2 (0,T;ffl(n))' 

where we used and ljT7)l . ■ 
Proposition 4.5 Let for any R £ 6 7?. £ emd ip,ip £ T>(Cl) 

^ = (l-tV Rc )ip + W RE G rE (^) (62) 

Then, for any r\ 6 2?([0, T[), we /iawe 

lim|$ e -^| n = (63) 

lim / p E u e § e rj{t)dxdt = / u<pr}' (t)dxdt + a vipr}' \t)dxdt. (64) 

Proof. The property l|63[) is an immediate consequence of (|35() and of the 
uniform boundness of G re (ijj) in L°°(f2). 
For the second property, let us notice that 

p e u £ <S> e -q'{t)dxdt = / / Xi\u e ^ 6 {x)rj'{t)dxdt 



n T Jo Jq 

+a E [ I u e G rs (ip)Tj'(t)dxdt. 
Jo Jd e 

As we obviously have 

lim / / Xfi u 6 § e (x)i]' (t)dxdt = / uiprf (t)dxdt, 
e ^°JoJn ' ' Jut 

it remains to study 

i-T i- i-T 



a. 



u e G rc (^)r]'(t)dxdt = a e \D £ \ f u E G re (ip)ri' \t)dxdt. 
o Jd e Jo J d e 



Using (|59|l and the uniform continuity of ip, we get 



£^0 



lim a £ / / u e Gr e (ij}) r q' '(t)dxdt = a / vtpr/'dxdt 



Proposition 4.6 If $> e is defined like in Proposition \4.5\ then we have 

lim f f Vu £ -V§ £ rj(t) dxdt = [ Vu-Vipr)(t) dxdt+4w^ [ (v-u)(tp-ip)rj(t) dxdt 
Jo Jn Jqt Jqt 

(65) 



Proof. First consider 

rT 



Vw e • V$ e dxdt 



o Jn, 
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which reduces to 

nVu £ ■ V^ry dxdt + [ [ Vu e • V$ e 7]dxdt. 
-h\c E Jo Jc E 

Lebesgue's dominated convergence theorem yields Vipln e \c s —> V</5 in L 2 (S1). 
Thus, taking (|4*§|l into account 

/ / Vw e • V(/??7Xa e \c e dxdt — > / S7u ■ S7(prj dxdt. 
Jo Jn Jn T 

Now, we come to the remaining part, namely 

/ / Vu £ -W^ e T](t) dxdt = f f (1- w Re )Wu e -Wipri dxdt 

JoJce Jo Jce 

I f (66) 
+ / / Vm e - Vw fle (G rt (t/0- ^)dxcft 

:= A+/ 2 

In the first integral, as Xc e V<^ — » in i 2 (Sl T ), Vw e — >■ Vw in L 2 {£l T ) and 
(1 — w Re ) is obviously bounded, we easily find that I\ tends to zero. 
In order to study I2, let us notice that 

I2 = I I Vu E • V% (G r (ip) — <p)r) dxdt+ 

Jo Jce ' (67) 

+ / / X7u e ■ X7w Rc {G r , O) - Gr E ((f)) V dxdt 

Jo Jce 

The first term in the right-hand side of i|t)7|) may be estimated by 

nVu £ ■Vw Re (ip - G rc (p))v dxdt\ < \Vu e \ nT \Vw Re ri\ n T\tp - G re (<p)|z°°(c e )- 

(68) 

As {w Re ) is bounded in ^ 1 (SI) (see Proposition 13. 2|) . the right hand side of l|68[l 
tends to zero by Remark 13.81 

Going back to the second term in the right hand side of (|67Jl . we may write 

fT r 

Vu e ■ Vw fli (G rE (i(;) - G re ((p))r)(t) dxdt 



JCc 



kez E V 



ip da - f tp da] I d$ / sin dQ 
si J ss 



T du e 
dr 



C k (r E M E ) " 



kez 

Airr E R 



V f t/jda-f tpda) {u £ \\ x - ek \ =rE - u £ \\ x _ £kl=R Jr)(t)dtdai 



( J R e - r £ ) fcc7 , y s; = e 7 s j e J J Sl Jo 

(G re (u £ ) - G Re (u E ))(G rc (^) - GrMHt) dxdt 



e 3 (R E - r e) j nT 
from which we infer that I2 is converging to 



4-7T7 / (v — u){i]-> — (p)ri(t) dxdt 
Jn T 
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and the proof is completed. ■ 

We are in the position to state our main result: 

Theorem 4.7 Thelimitsue L°°(0, T; L 2 (Q))nL 2 (0, T; ^(fl)) andv E L 2 (Q T ) 
°f \4tf) verify (in a weak sense) the following problem: 

^ - Ait + 4tt7(u -v) = (/ - g) in fF, (69) 

a— + 4:irj(v — u) — g in f2 T , (70) 

«(0) = wo in Q. (71) 
u(0) = v in CI (72) 

Moreover, there holds u £ C°([0, T];L 2 (fl)) and v € C°([0, T]; i/" _1 (Q)); i/iese 
are i/ie senses of PHj) and J7ij| ). 

Remark 4.8 ^4s i/ie problem i)6'^) - |7S| ) /ias a unique weak solution, the conver- 
gences in Proposition ^. g| feo/rf on ifte whole sequence. 

Proof of Theorem 14.71 We set in (|17|l u> = $ e where <£> e is defined like in 
lemma l4~51 Then, by multiplying (|17fl by 77 <E T)([Q,T[) and integrating it over 
[0, T] we get 

p £ u £ $ £ r 1 'dxdt + [ k £ \7u £ (\7$> £ ) 71 dxdt = [ (f £ ,$ £ )ndt + [ p £ u £ <S> e r]{0)dx. 
Jq. t Jo Jn 

(73) 

Then, the left-hand side tends to 

— / utprj'dxdt — a, vipn'dxdt + / Vu ■ Wip r\ dxdt+ 
Jn T Jn T Jn T 

(74) 

+47T7 / (v — u)(ip — ip) n dxdt. 
This is a direct consequence of Proposition 14. 61 together with the remark that 

[ [ Vu £ V$ £ dxdt = 

JO J D e 

since $ e is constant on every B(ek, r £ ), k € Z £ . 
As for the right-hand side, we have 

(f £ 7 <Z> £ ) v dt = f (F,(l-w Re )tp)r ) dt+ [ (f £ ,w RE G r M))vdt 



T 



and, with hypothesis i|46[l . 



(f £ ,(l~w R J<p}r,dt^ / (f-g,<p)v& 



Moreover, 

J (f £ ,w Rc G rc ^)) V dt= j (/••.«'„■ ((V,. ( < -)-r.)/,r//+ / (/••.«•,,. r),,<ll 
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with 



T (f,w Rs (G r M-^))vdt 







• T 

< I \f e \ H -^R £ (Gr e W - n H nny 



As we have 

\wR.{G r .{1>) - $)\ttkm = \V[w R AGrM - ^))|n 

< \Vw Re \ Ce \G rc (V) - ^U«.(C.) + I V^kuD, 

Remark 13.81 and l|34|) obviously yield 

lim|w,R £ (GV e (V>) - VOI^n) = 0. 
The assumption (|2U|) on / 6 implies that |/ e |_H--i < C and thus 

Urn / (f,w Re {G r M - il>))t]dt = 0. 



We conclude thanks to hypothesis (|46|l that 

lim / (f e ,w Re ip)r]dt = / (g, tp)r]dt. 
£ ^°Jo Jo 

Finally: 

lim / (f,$ s )r 1 dt=[ (f-g,<p)r 1 dt+ [ {g^)r)dt. 
E ^°Jo Jo Jo 

We get 

/ p £ u e <f> E r)(0)dx = / u E o <S> e n{0)dx + a e u E G rE (^)r](0)dx. 
Jn Jn E J d e 

Using the hypotheses (|21() (|23H on Uq, we pass to the limit and with the same 
arguments as above we obtain 

lim / p e u e ) & e r)(0)dx = 7/(0) / (u ip + avoip)dx 
Jn Jn 

which achieves the proof. ■ 

5 Homogenization in the case e 3 << r £ « e 

In this section, we fix some R £ € lZ e . 

Remark 5.1 Notice that in this case Proposition \3. 7| also reads 

\tp\ 2 dxdt < C\Vcp\ 2 L2(nT) , Vcp G L 2 (0,T;Hl(Sl)). (75) 



1 

'0 J D E 



In the present case, Proposition 12 . 31 and Lemma |3. til imply in a straightfor- 
ward manner the result corresponding to Proposition ^. 21 



1G 



Proposition 5.2 There exists u£ L°°(0, T; L 2 (Q))nL 2 (0, T; Hq(Q)) such that, 
on some subsequence, 

u e ^u in L°°(0,T;L 2 (Q)) (76) 
u e ^u in L 2 {Q,T;Hl{n)) (77) 
G fle (u £ )-*« in L 2 (n T ) (78) 
G re (u £ )^u in L 2 (fi T ) (79) 

Moreover, we have 

lim / / \u s -G r (u e )\ 2 dxdt = (80) 
Jo J D c 

In the light of Remark 15. II we prove as in the previous section: 
Proposition 5.3 For any tp £ L 2 (0, T; Hq (Q)) ; we have 

/ +■ u e ipdxdt — > / utfdxdt. (81) 
Jo 7 D e 7n r 

The homogenization result obtained in this case follows. 

Theorem 5.4 The limit u € L°°(0,T;L 2 (n))r\L 2 (0,T;H^(n)) of JTZty-fT^ is 

the only solution of 

(l + o)^-Au = / in ft T , (82) 

u{0) = jrhrf + jiTaj v " m n (83) 

Moreover, the convergences in Proposition 15.21 ftoZd on </ie whole sequence and 
u £ C°([0,T];Z 2 (fJ)), this being the sense of f%fy . 

Proof. The proof of l|82|) is similar to the corresponding one of the Theorem l4.7l 
The test function <f> £ is given by 

$ e = (i -w Re )<p + w RE G rc (<p), peP(ft). 

The only interesting convergences are the following two: 
\7u £ (\7wR c )(G rE (ip) — ip)dxdt 

CT 



< C|Vw e | T|Vw fl J T \G re ((p) ~ V| L ~(CT) ^ 

<c 7 y^.c(^) 1/2 (f)^o 



rT {f s ,w R AG r M-v)) 



< C\(G re (<p) - (p)Vw Re \ L 2 {c t } +C\w Re X7<p\ l2{CTuDT) 



< 
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< c \v<p\ LX(u) ( 7s 1/2 i4 + \e s u d^ 1 ' 2 ) - o, 

where we have used the a priori estimates of Proposition 12.31 Remark 13.81 and 
Proposition ^. 21 

Using Proposition !,^ 21 and hypotheses H21|) - (|23[l . we obtain with the same argu- 
ment as before 

lim / p e u e a Q e n(Q)dx = 77(0) / (u + av )ipdx 
Jn Jn 

which achieves the proof. ■ 



6 Homogenization in the case r £ « e 3 . 

As in this case 7 e — > 0, we only can prove: 

Theorem 6.1 There exists u £ L°°(0, T; L 2 (n)) n L 2 (0, T; H&(Sl)) such that 

u £ ^u in L°°(0,T;L 2 {n)) (84) 
u £ ^u in L 2 {0,T;H^(Q)) (85) 

u e X D e ^v in a.e. fe[0,T] (86) 



\D L , 

where u is the only solution of the following problem: 

^-Au = f in Q T (87) 
at 

u(0) = u a in il (88) 

Proof. The convergences (|84|) - (|85|l hold on some subsequences; they are in- 
sured by Proposition 12. 31 We have to remark that (|25|l is the hypothesis which 
insures the existence of v G L°°(0, T; L 2 (£l)) which satishes 

-^u £ X D E ^v in V\n) a.e. te[0,T] 

(see Lemma A-2 j^]); we have to prove that v = vq. 
Acting as usual, we take 

$ £ - (l-w Rs )ip + w Re G r M (89) 

for some R £ £ lZ e and ip, i/j £ 2?(S1). Notice that in this case we have 

<5> e ^ip in H%(Q) (90) 

because obviously wr s —>■ in Hq(Q), 

Passing to the limit in the variational formulation, we obtain in a straightforward 
manner 

— / wpn' dxdt — a I viprj l dxdt+ / \7u\7ipi]dxdt = / (f,Lp)r/dt+ 
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+ ( / u <fdx + a f vaipdx) T)(0), V77 S X>([0,T[) 

Setting if — 0, we find that v is independent of t and that v £ C°([0, T]; L 2 (Q)), 
which achieves v = vq. Then, setting ip = 0, we prove (|57jl and l|5S|l. the last 
one holding also in the sense of C°([0, T]; L 2 (il)). ■ 
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